Addressing atoms in optical lattices with Bessel beams 
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A method of synthesizing localized optical fields with zeroes on a periodic lattice is analyzed. The 
applicability to addressing atoms trapped in optical lattices with low crosstalk is discussed. 
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There is much current interest in using atoms 
trapped in optical lattices for quantum logic 
devices 0,0. One of the challenges in implementing 
this scheme is posed by the need to address indi- 
vidual atoms with near resonant light. Current esti- 
mates of error limits for scalable quantum computing 
require primitive logic operations with errors as low 
asjE Q O(10- 5 - 1(T 6 ). The corresponding limita- 
tion on optical crosstalk when addressing one atom 
in a lattice can only be quantified in the context 
of a particular choice of quantum gate. In atomic 
schemes error rates for operations that are depen- 
dent on single photon processes tend to scale with 
the intensity. This implies that the intensity leakage 
of a logical control or state readout beam at a site 
adjacent to a site being addressed should not exceed 
0(1O -5 ). In this letter I describe a novel approach 
to image formation that allows sites in an optical 
lattice to be addressed with minimal crosstalk. 

In order to emphasize the difficulty of achieving 
low crosstalk in optical lattices consider the follow- 
ing simple scaling relationships. We have a 1- or 
2-D lattice lying in the x — y plane that is defined 
by countcrpropagating beams at wavelength A / . In- 
dividual atoms are separated by a minimum dis- 
tance of d — A//2 and we wish to address them 
with near-resonant light of wavelength A. The most 
obvious approach to doing so involves focusing a 
beam propagating perpendicular to the plane of 
the lattice (along z) to a small spot. Assuming 
a Gaussian beam profile the intensity distribution 
is I{p,z) — I exp(— 2p 2 /w 2 ) where p 2 = x 2 + y 2 , 
w 2 (z) — w 2 (l + z 2 /z R ), wq is the beam waist at 
z = 0, and zr — ttwq/X. Expressing the ratio of the 
intensity at a neighboring site to the on-site inten- 
sity as e implies a beam waist wo = \f— l/(21ne) A/. 
For e = 10~ 5 this evaluates to w = Wo/Xf = 0.21. 

We can quantify the corresponding requirement 
on the ratio of trapping and addressing beam wave- 
lengths in terms of the performance of high numer- 
ical aperture lenses. We assume a lens system with 
aperture diameter D is used to focus the addressing 
light. If the ratio of the lens aperture to the Gaussian 
waist at the aperture is given by p = D / w{z\ Dns ) then 
a focal plane waist of wq implies D ~ pwoZi ens / zr. 




FIG. 1: Lens numerical aperture needed to focus a Gaus- 
sian beam to a waist wq for A//A = 1, 2, and 10. 



A reasonable minimum requirement on p to avoid 
rings surrounding the focal plane image is that the 
fraction of the optical power blocked by the aper- 
ture should not exceed the desired intensity ratio e 
by more than a few times. Using an optimistic value 
of p = 3 gives approximately 1% transmission loss 
at the lens and a numerical aperture of 
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The variation of the lens NA with wq is shown in Fig. 
^ Available microscope objectives, as well as high 
resolution optical lithography lens systems, tend to 
have numerical apertures not more than ~ 0.8 un- 
less oil immersion is used which is not compatible 
with atomic imaging inside vacuum chambers. This 
implies that the lattice light wavelength A / must be 
several times longer than the addressing wavelength 
A in order to address single atoms with negligible 
leakage to neighboring sites. Thus a lattice of Rb 
atoms which can be manipulated with light near res- 
onant with the D2 line at A = 0.78 /xm will not 
be individually resolvable in the context of scalable 
quantum logic unless A/ > 1.5p.m. 

Various solutions to the addressability problem 
are under study. The most direct solution is to 
make A/ 3> A. Arrays of widely spaced traps 
that use many diffractively generated beams in- 
stead of lattices also fall into this category of solu- 
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tion. Unfortunately this is not compatible with load- 
ing from a Bose Einstein condensate via the Mott 
transition^ , although other loading schemes may 
still be usedQ, ||- Another possibility is to keep 
A ~ Xf but load the lattice so that only every few 
lattice sites are occupied, or change the angle be- 
tween the lattice beams after loading so that a longer 
periodicity is obtained^. Hc|. 

Here we analyze an alternative approach to sin- 
gle atom addressing in optical lattices with A ~ A/. 
The idea is to accept the fact that the light cannot be 
localized sufficiently well, but tailor the beam pro- 
file so that the field has zeroes at neighboring lat- 
tice sites. The basic geometry is sketched in Fig. 
121 N beams, each with a phase controlled by a 
spatial light modulator pixel, propagate in the xy 
plane and are aligned to converge on the origin at 
x — y — 0. The beams necessary for creating the 
optical lattice can in principle be combined with the 
imaging beams using dichroic mirrors. Each con- 
verging plane wave is polarized along z with am- 
plitude Aj exp[i(— kj • p — tot + Xj)} + c.c., where 
k 3 = k(cos<f>jX + smcj)jy), k = 2n/X, <pj is the az- 
imuthal angle of wave j, and Aj, Xj are adjustable 
amplitudes and phases. As shown in the figure this 
can be achieved using a single one-dimensional mod- 
ulator with N pixels and a system of mirrors and 
lenses. In the limit when N — > oo and the beam am- 
plitudes and phases are all equal the field amplitude 
generated on axis is just the zero order Bessel beam 
J (kp). 

When the field is synthesized from a finite num- 
ber N of plane waves it is quasiperiodic in space, and 
rings of secondary intrerference maxima occur as can 
be seen in Fig. [21 The diameter of the secondary 
rings can be estimated from Afcd r i ng ~ 2n. Using 
Afc = 2k /(N/2) gives d ring ~ NA/4, while the ac- 
tual diameter found numerically is some 25% larger. 
Recalling the shift theorem of Fourier analysis it is 
straightforward to scan the spot location over a dis- 
tance of up to d r ing by adding appropriate phase off- 
sets to the incident beams, as is shown in Fig. |21 It 
is possible in this way to address hundreds of atoms 
using a multipixel modulator, without mechanical 
motion of the optical system. This approach to syn- 
thesis and scanning of localized optical fields also 
finds application in atomic lithography [ll|. 

A Bessel beam written in this way with A = 
0.78 /im has a central lobe with 1/e 2 intensity ra- 
dius of 0.22 pm. However the closely spaced sec- 
ondary maxima of the Bessel function give crosstalk 
on nearby lattice sites. We describe now a method of 
synthesizing a beam with zeroes at regularly spaced 



sites of a one-dimensional lattice. An arbitrary solu- 
tion of the two-dimensional Hclmholtz equation can 




FIG. 2: Optical layout for writing a Bessel beam. The 
intensity distribution was calculated with N = 100 and 
A = 0.78 /Urn. The central lobe was translated by Aa; = 
4 pm and Ay = 2 pm by adding phase offsets, see the 
text for details. 



be written as a Fourier-Bessel series 

oo 

A(p, 9) = J (fcp) + a n J n {kp) exp[in6] (2) 

71=1 

with the a n complex coefficients to be determined. 
We have set the zero order coefficient to unity corre- 
sponding to a normalized field with unit amplitude 
at the origin. Any field of the form J2J) can be syn- 
thesized with an amplitude and phase modulator in 
the geometry of Fig. [21 To make this explicit we 
note that the Fourier transform of (J21 is 

/>oo t*2ir 

A(q,<t>) = I dp p d9A(p, 6) exp[ipqcos(9 — 4>)] 
Jo Jo 

= —S(q - k) (l + a ni n exp[m$ J . (3) 

Thus a ring of converging plane waves with 
wavenumber k and complex amplitudes given by the 
term in parentheses in Q provides the desired field. 
The role of the lenses shown in Fig. [21 is to tilt the 
wave passing through each spatial light modulator 
pixel towards the center of the image region, but 
the lenses do not tightly focus each wave. Outside 
of the central spot destructive interference leads to a 
low background light level. From the point of view 
of Fourier optics focusing of a beam results from 
interference of the constituent plane wave compo- 
nents, taking into account relative phase shifts due 
to propagation. The usefulness of the approach to 
image synthesis presented here is that we can di- 
rectly control the amplitudes and phases of each of 
the plane wave components. 
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TABLE I: Bessel coefficients and resulting maximum crosstalk which occurs at site m max away from the origin for 
A = 0.78 /jm and A/ = 0.8 /im. The last row shows the maximum error in the first 50 neighboring sites using expansion 
© with N = 256. 



A field consisting of a finite number of terms in 
© that is useful for addressing a one-dimensional 
lattice along x will have A(p, 0) = A(p, ir) which im- 
plies that a n = for n odd. Requiring that A vanish 
at the M lattice points specified by p m = mXf/2, 
6 m = 0, for rn — 1...M and limiting the sum in J5J 
to M terms gives M linearly independent equations 
Jo(kpm) + Sn=i a>2nJ2n{kp m ) = 0. These equations 
are easily solved for the coefficients ci2n- The coef- 
ficients decrease rapidly with Bessel order provided 
Xf is not too much less than A. The series coeffi- 
cients, and maximum intensity crosstalk at any site 
not being addressed are listed in Table |U We see 
that M = 6 is sufficient to ensure a crosstalk of a 
few times 10 -5 . In an actual implementation there 
are several limiting factors to consider including the 
number of spatial light modulator pixels, and the 
amplitude and phase resolution of each one. The 
last row in Table [I] shows the crosstalk is up to sev- 



eral times higher than the theoretical value if we as- 
sume the field is synthesized from 256 beams equally 
spaced azimuthally, with 14 bit resolution in ampli- 
tude and phase modulation for each beam. Addi- 
tional calculations with A/ = 1.0 (im give crosstalk 
levels roughly 10 times lower than the example in 
Table IJ 

In summary we have described a novel method of 
addressing atoms in periodic 1-D lattices with low 
crosstalk. The extension of this approach to a 2-D 
lattice is complicated by the fact that a 2-D lat- 
tice has angle dependent interatomic spacings. The 
Fourier-Bessel expansion (J2J) can still be used, but 
the coefficients tend to grow rapidly with n. Gener- 
alizations that are suitable for a 2-D geometry are 
currently under investigation. 
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